Sound speed of a Bose-Einstein condensate in an optical lattice 
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The speed of sound of a Bose-Einstein condensate in an optical lattice is studied both analytically 
and numerically in all three dimensions. Our investigation shows that the sound speed depends 
strongly on the strength of the lattice. In the one-dimensional case, the speed of sound falls mono- 
tonically with increasing lattice strength. The dependence on lattice strength becomes much richer 
in two and three dimensions. In the two-dimensional case, when the interaction is weak, the sound 
speed first increases then decreases as the lattice strength increases. For the three dimensional 
lattice, the sound speed can even oscillate with the lattice strength. These rich behaviors can be 
understood in terms of compressibility and effective mass. Our analytical results at the limit of weak 
lattices also offer an interesting perspective to the understanding: they show the lattice component 
perpendicular to the sound propagation increases the sound speed while the lattice components 
parallel to the propagation decreases the sound speed. The various dependence of the sound speed 
on the lattice strength is the result of this competition. 



PACS numbers: 03.75.Fi, 03.75.Kk, 05.30.Jp 

I. INTRODUCTION 

A Bose-Einstein condensate(BEC) in an optical lattice 
has attracted great interests recently, both experimen- 
tally and theoretically [U, [1]. The presence of the lattice 
can remarkably enrich the behaviors of the system com- 
pared to the uniform case, providing a new fertile ground 
for exploring a variety of solid-state effects in BECs, for 
example, Bloch oscillations [1, 0, Q and Landau-Zener 
tunneling!^, between Bloch bands in an accelerat- 

ing optical lattice. Moreover, a BEC in an optical lattice 
can be considered as quantum simulators and therefore 
used for testing fundamental theoretical concepts 01 ■ For 
example, it can be used to simulate the Bose-Hubbard 
model and study experimentally the quantum phase tran- 
sition between superfluid and Mott-insulatorp^ [Tlj . 

In this article, we launch a systematic study of the 
speed of sound for a BEC in an optical lattice in all three 
dimensions. The speed of sound is important for two 
simple reasons: (1) it is a basic physical parameter that 
tells how fast the sound propagates in the system; (2) it 
is intimately related to superfluidity according to Lan- 
dau's theory of superfluid. Because of these, the sound 
propagation and its speed was one of the first things that 



have been studied by experimentalists on a BEC since its 
first realization in 1995 The propagation of sound 
in a harmonically trapped condensate without the lat- 
tice has already been observed experimentally [isj and 
studied theoretically [14]. Now there are experimental ef- 
forts to measure the sound speed for a BEC in an optical 
lattice [H. 

There have been a great deal of theoretical work done 
to understand the sound speed for a BEC in an opti- 
cal lattice. These studies show that three parameters 
strongly affect the speed of sound: the strength of the op- 
tical lattice V, the interaction between atoms c, and the 
lattice dimension D (D = 1, 2, 3). In Ref.fT^, the phonon 
excitations of the BECs in one-dimensional optical lat- 
tice [D = 1) are theoretically investigated by solving 
the Bogoliubov equations. Their analytic results for the 
sound speed in the weak potential limit predicted that 
the sound speed decreases monotonically with increas- 
ing the depth of the optical lattice. The most detailed 
study of sound propagation in one-dimensional lattice 
was done by Pitaevskii and Stringari's group [l7l[l8|. who 
also found the sound speed is suppressed by the lattice. 
In particular, Ref. presents the detailed comparison 
between the sound speed obtained by Bogoliubov theory 
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with the one obtained from the compressibihty and the 
effective mass. Similar resultsfl^ were also obtained for 
the Kronig-Penney potential, a special form of the peri- 
odic potential. Furthermore, Martikainen and Stoof'2^ 
examined the effect of the transverse breathing mode on 
the longitudinal sound propagation for a BEC in a one- 
dimensional optical lattice. In particular, they discuss 
how the coupling with the transverse breathing mode in- 
fluences the sound velocity in an optical lattice. Kramer 
et al. [2l| also studied the effect of the transverse degrees 
of freedom on the velocity of sound of a BEC in ID op- 
tical lattice and radially confined by a harmonic trap. 
A recent paper by Taylor and Zaremba[22| studied the 
Bogoliubov excitations of a BEC in an optical lattice in 
all spatial dimensions. However, the formulation in Rcf. 
[221] is in principle; the authors did not present the con- 
crete results of sound speed in two and three-dimensional 
cases(D=2,3). Most interestingly, with numerical calcu- 
lations Boers et al.^^ found that the sound speed of 
a BEC in a three-dimensional optical lattice achieves a 
maximum with increasing lattice depth. Because of the 
difficulty to obtain the Bloch states with interaction, the 
investigation of Boers et al. is limited to the low density 
so that the Bloch wave function of the free particle can 
be used as an approximation. 

Our investigation here tries to overcome the deficien- 
cies in previous studies to give a complete picture how 
the sound speed is affected by the lattice strength u, 
the interaction between atoms c, and the dimensional- 
ity D. Analytical approaches are used in two limiting 
cases: weak lattice and strong lattice. For weak lattices, 
they can be viewed as perturbations. In this case, we ob- 
tain an analytical expression to the second order of the 
lattice strength for the sound speed of a BEC in an ar- 
bitrary periodic potential. We have analyzed this result 
for the important case of the periodic potential being an 
optical lattice. Our analysis finds a strong dependence 
of the sound speed on the lattice dimensions. Especially, 
we find that the lattice component perpendicular to the 
sound propagation increases the sound speed while the 
lattice components parallel to the propagation suppresses 
the sound speed. Since the lattice can only be parallel 
to the propagation direction of sound in one dimension 
{D = 1), the sound speed falls monotonically with in- 
creasing lattice strength. In two and three dimensions 
{D = 2,3), there are both perpendicular and parallel 
components in the lattice and, therefore, there is a com- 
petition. As a result, there is a rich dependence of the 
sound speed on lattice strength in the case of D = 2,3. 
The sound speed can first increases then decreases as the 
lattice strength increases. We have also tried to under- 
stand these results from a different angle, i.e., in terms 
of compressibility k and effective mass m* . The analyti- 
cal expression is found for compressibility k and effective 
mass m* for a BEC in an optical lattice. We find that the 
effect of the lattice on the sound speed reflects a competi- 
tion between the slowly decreasing compressibility k, and 
the increasing effective mass m*, with increasing lattice 



depth. 

In the limit of strong lattices, it is reasonable to use 
the tight-binding model to describe the BEC in an optical 
lattice^]- Our analytical results display that the sound 
speed always exponentially decreases with increasing the 
optical lattice in all dimensions. This D-independent be- 
havior of sound speed can be understood as the competi- 
tion between the tunneling strength J between adjacent 
sites and the interaction U between the atoms at a lattice 
site. With increasing lattice depth, U slowly increases 
while J exponentially decreases, resulting in monotoni- 
cally decreasing speed of sound. 

Our analytical results are complemented by our numer- 
ical study, where the results are obtained for all ranges 
of lattice strength. Our numerical results agree well with 
our analytical results in both weak potential and tight- 
binding limits for the case of weak interatomic interac- 
tion. For the intermediate strength of lattice, we find that 
the sound speed even oscillates with the lattice strength 
for a three-dimensional optical lattice. We emphasize 
that in our numerical calculations the interaction be- 
tween atoms is taken into account to compute the Bloch 
states in all three dimensions. In Ref.jl^, the interaction 
is neglected in computing Bloch states for BECs and the 
Bloch states of free bosons were used as an approxima- 
tion. 

This paper is organized as follows. In Section II, for 
the sake of self-containment and introducing notations, 
we describe the basic theoretical framework of our study. 
It includes the definition of the sound speed Vs^ com- 
pressibility K, and effective mass m* . In Section III, we 
present the analytic results of the sound speed for a BEC 
in the optical lattice in both weak potential limit and 
tight-banding regime. Section IV contains our numerical 
study of the sound speed. The details of our numerical 
methods are given here. In Section V, we discuss the 
possibility of observing the phenomena presented in this 
paper within the current experimental capability. The 
last section (Sec. VI) contains a discussion of our re- 
sults and concluding remarks. Five appendices are given 
at the end to show the detailed steps to derive our key 
analytical results in the main text. 



II. BASIC THEORY 

A. Mean-field theory of Bose-Einstein condensates 

We focus on the situation that the BEC system can be 
well described by the mean-field theory. In this case, the 
BEC system is governed by the following grand-canonical 
Hamiltonian, 



d^r 



{r (r) 



1. 



Vlatt{l^\4'{r) 

|^(r)|^}. (1) 
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In our case, the external potential is a three-dimensional 
optical lattice created by six laser beams that are per- 
pendicular to each other]!], [3| , 



Via 



(f) = cos {x) -I- cos (y) -I- cos (z) 



(2) 



where v characterized the strength of the optical lattice. 
In Eq. ([T]), all the variables are scaled to be dimensionless 
by the system's basic parameters: the atomic mass m, 
the wave number oi the laser light, and the average 
density tiq. The chemical potential ^ and the strength v 
of the periodic potential are in the units of Ah^k^/m, the 
wave function ij) is in the units of y^no, r is in the units 
of l/2fci, and t is the units of m/Ahk\. The nonlinear 
coefficient c — T:n^as/k\, where > is the s-scattering 
length. 

Sound is a propagation of small density fluctuations 
inside a system. To study sound in a BEC, one first 
need to find out the ground state of this BEC system, 
which serves as a media for sound propagation. The 
sound speed can then be found by perturbing the ground 
state as explained in detail in the next subsection. 

The ground state of a BEC in an optical lattice is a 
Bloch state at the center of the Brillouin zone. Briefly, 
the Bloch state is of the following form 



(3) 



where k is the Bloch wave vector and <f>j: {r) is a periodic 
function with the same periodicity of the optical lattice. 
The Bloch wave function 4>^{r) satisfies the following sta- 
tionary Gross-Pitaevskii equation 



\{y + ^k)' 



where is the chemical potential. The energy of the 
system in a Bloch state is given by 



E{k\ = 



W + ik 



The set of energies E{k) then forms a Bloch band|25l.[2q. 
The Bloch state can be obtained analytically in certain 
circumstances|27l. In most cases, it has to be computed 
numerically [25, 28]. The numerical method of this study 
is described in Section VI. To compute the sound speed, 
one may only need the Bloch state (po a.t k = 0. However, 
for the effective mass defined by 



1 

m* 



dk^ 



(6) 



where k = \k\, one also has to compute Bloch states in 
the vicinity of fc = 0. In this article, we also study one- 



and two-dimensional cases. The one-dimensional optical 
lattice is given by 

V{x) = ucos(a;) , (7) 
and the two-dimensional optical lattice is given by 

V{x, y) — v[cos{x) + cos(j/)] . (8) 

B. Definitions of the sound speed 

In Section II. A, the BEC system is regarded as a 
Hamiltonian system by the grand canonical Hamilto- 
nian ([1]); the corresponding Gross-Pitaevskii equation 
can be obtained by the variation of the Hamiltonian, 

dtb 1 

i^ = --^^^ + V{r)^P + c\i,\^ij. (9) 

The Bogoliubov equations can be determined from the 
linear stability analysis of the CP equation ([9]). To ex- 
plore a small disturbance S(f)j^ (r , t) at a Bloch state (pj: (r ) , 
we write 

V,(r, t) = e^-^-^^^' [(/.J (f) + % (f, t) ] , (10) 
where the disturbance can be similarly written as 

SPj: (r, t) = u^e*{«-*'-'('?^*> + y*_.e-i{?-^-<9)t} ^ (n) 

Plugging Eq. (fTTjl into Eq. ([9]) and keeping only the 
linear terms, we arrive at the Bogoliubov equations [25l|, 



with 




^{k + q 



and 



1 

-1 



(12) 



(13) 



(14) 



where ^ {q) is defined as 

^ (q) - -^(V + iq)' + V{r)-^l + 2c|0,,|2. (15) 

Note that q represents mode of the small perturbations 
and is of nature of Bloch wave vector as the matrix M is 
periodic. 

In general there are two equivalent definitions for the 
sound speed in a BEC. As sound can be regarded as a 
long wavelength response of a system to a perturbation, 
the sound speed can be extracted from the excitation of a 
BEC. According to the Bogoliubov theory, the excitation 
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energy e {q) of the BEC in a Bloch state at fc = can be 
found by solving the eigenvalue problem of Eq. (H^)- In 
terms of these excitations, the sound speed of a BEC 
system can be defined as 



(16) 



where g = |(7|. 

The other definition arises when the BEC system is 
regarded as a hydrodynamics system. In this context, 
the sound speed in a BEC is given by the standard 
expression[l3, [H ^ 



1 



(17) 



where k is the compressibility of the BEC system and 
defined as 



no 



dno' 



(18) 



where the chemical potential fi and no is the averaged 
density. For a BEC system with repulsive interatomic 
interaction, the optical trapping reduces the compress- 
ibility of the system as the effect of the repulsion is en- 
hanced by squeezing the condensate in each well. Ac- 
cording to the definition of sound speed in Eq. p7|) , the 
sound speed refiects the competition between the com- 
pressibility K and the effective mass m* . 

Both definitions are used in our computations and they 
agree with each other as expected. The proof of the 
equivalence of these definitions can be found in Ref. [2^ . 



III. ANALYTICAL RESULTS 

A. Weak potential limit. 

We consider first an arbitrary periodic potential Var (f) 
with the periodicity of R, 



with 



R = niiai + 771202 + m^as, 



(19) 



(20) 



where r is the position vector, di, 02, and 03 are the three 
primitive vectors and mi, m2, 7713 range through all in- 
tegral values. In the weak potential limit, the periodic 
potential Var (r) can be regarded as a perturbation. This 
allows us to solve both the Gross-Pitaevskii equation (|3|) 
and the Bogoliubov eigenvalue problem (jl2p perturba- 
tively by expanding the wave function and chemical 
potential /i of the BEC system in the order of the weak 
potential. 



7^(1) 
(1) 



4- „(2) i 



where is zeroth order of the potential strength, 

/i'-^-' first order, etc. We find that the sound velocity 
along a given direction indicated by a unit vector f is 



n 



,}^n'{V). (22) 



\n ■ 



Inl^iAc+lnl^y 



In the above, ^n{V) is the Fourier coefficient of Var {"r) 
as defined by 



with 



■^n{V) = i / dWa.(f)e- 
" J cell 



n = nibi + 712^2 + n^b^ 



(23) 



(24) 



where nj^s are integers and &j's are the set of reciprocal 
primitive vectors defined by 



Gi ■ bj = 271 Si- 



(25) 



In the integration, is the volume of the primitive cell 
and the integration is over one primitive cell. The de- 
tailed derivation of Eq. can be found in the Ap- 
pendix C. 

The focus of this article is optical lattices as described 
in Eas. (|2l7l8p . In this special but important case, the 
primitive vectors oi, 02, and 03 can be chosen along the 
directions of laser beams, x, y, and z, respectively. Also 
we have |ai| = |a2| — l^sl = Stt. For this case, we find 
from Eg. ip^ that if the sound propagation direction is 
along the x-axis, the sound speed is (see also Eq. (jC10[) 
in the Appendix C or Eq. (|D43p in the Appendix D) 



V. = Vc - 



E 



81/c {nl + nl) |np - Acn\ 



|n|2(4c-f |ri|2) 



2^3 



^r\V).{2Q) 



(21) 



The sound speeds along the y-axis and z-axis can be 
found easily with permutation argument and the sound 
speed along a general direction is a certain combination 
of these three speeds. 

When there is no periodic potential Var (r) — 0, the 
sound speed in Eq. (|26p is reduced into ^/c, the sound 
speed for a BEC in free space, as expected. We also no- 
tice that there is no first-order correction to the sound 
speed due to the periodic potential. Most importantly, 
the analytical result in Eq. reveals that the lat- 

tice component perpendicular to the sound propagation 
(generated by the laser beams along the y and z-axes) 
increases the sound speed while the lattice components 
parallel to the propagation (generated by the laser beams 
along the x-axis) decreases the sound speed. As a result 
of this competition, the sound speed can either increase 
or decrease with lattice strength. This competition be- 
tween the parallel and perpendicular components of the 
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optical lattice certainly also applies to a general periodic 
potential if one carefully examines Eq. (j22p and interprets 
"parallel" and "perpendicular" in a more sense. 

To illustrate this more clearly, we consider a sim- 
ple case where the periodic potential is a ID optical 
lattice given by Vexir) = vcos{y). There are only 
two non-vanishing Fourier coefficients, i.e. ^0,1. o(^) = 
= vf^- Then according to Eq. the 



speeds of sound along the x, 
tively. 



y and 2-axis read, respec- 



1 



f - 



(27) 



which show that with increasing the strength of the op- 
tical lattice the sound speed along the y-axis, parallel to 
the periodic lattice falls while the speeds of sound along 
both the X and z-axes increase. 

Now we study the BEC sound speed in optical lattices 
in terms of compressibility and effective mass according 
to the second definition of speed of sound, i.e. Eg. lfTT)) . 
Again we treat the weak optical lattice as a perturba- 
tion. For optical lattices of all dimensions as described 
in Eqs. (|2l7l8p . we find the chemical potential at fc = 



4 (2c 

and the system energy near k ~ 



(1 + Acy (4fc2 - 1) 



(28) 



(29) 



So, the chemical potential depends on D, the dimension 
of the lattice, while the energy E{k) does not. The com- 
pressibility K can be calculated from the chemical poten- 
tial iJL via Eq.p^ and it is given by 



Dcv^ 



2c)- 



(30) 



This shows that the compressibility k tends to decrease 
with increasing v as the optical lattice localizes the BECs 
inside each well. Moreover, the compressibility k de- 
creases faster with v for higher dimensional lattices. The 
effective mass m* can be computed from E{k) and it is 
found that 



1 

m* 



= 1- 



2v' 



2c 



(31) 



dependence of the effective mass on v is independent of 
the lattice dimension D. As the speed of sound is de- 
fined as Vs = 1 / (nm*) , the compressibility k and the 
effective mass m* influence the sound speed in opposite 
directions. 

Plugging both Eqs.jSOl) and jSl]) into Eq. (HH), we find 
an analytical expressions for the sound speed of a BEC 
in an optical lattice up to the second order of v 



V (4c +1)' , 



(32) 



With simple algebra (see Eq. (jE6|) in Appendix E), one 
can show that this expression is consistent with the more 
general formula in Eq.f^H). 

Eq. p2| indicates that in one dimension {D — 1) 
the effective mass m* always wins over the compress- 
ibility K, resulting in decreasing sound speed with the 
lattice strength. However, in two or three dimensions 
{D — 2,3), the situation is very different. There exists 
a critical value of c, the interatomic interaction strength, 
beyond which the effective mass m* wins. Otherwise, 
the compressibility k has bigger influence on the speed 
of sound and the speed of sound increases as the lattice 
becomes stronger. The critical values are c = 1/4 for 
D = 2 and c = 1/2 for D = 3. 

We have discussed the behavior of the speed of sound 
in two different languages: one in terms of perpendicular 
and parallel components of the periodic potential with 
Eg. ipS)) and the other in terms of effective mass m* and 
the compressibility k. Are these two pictures consistent? 
The answer is yes. To see this, we re- write Eq. ipS)) as 



(4c-f |n|2)3 



E 

n=iO 



n\^{4:C+ |n|2) 



^}^n\V). (33) 



By comparing it with Eqs. (|5D|) and (PT|) . it is apparent 
that the first term in the curly brackets comes from the 
compressibility k while the second term results from the 
effective mass m*. This observation gives us the follow- 
ing picture: on the one hand, all components of the peri- 
odic potential contribute to the compressibility k, which 
leads to the enhance of the sound speed; on the other 
hand, only component parallel to the sound propagation 
increases the effective mass m*, which leads to the sup- 
pression of the sound speed. Since the effective mass m* 
always wins over k along the parallel direction, we come 
to the previous understanding: the perpendicular com- 
ponents increase the sound speed while the parallel one 
suppresses it. 



It is clear that the effective mass always increases with 
the lattice strength v. This is expected as the increased 
lattice strength suppresses the tunneling between neigh- 
boring wells thus increases the effective mass m*. In- 
terestingly, in contrast to the chemical potential /i, the 



B. Tight-binding regime 

When the potential wells are sufhciently deep, the con- 
densate is well localized at each lattice site and the fol- 
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lowing tight-binding model may become adequate to de- 
scribe the BEC in an optical lattice [2^ 



<n,n'> n 



where the first summation is over all pairs of the nearest 
neighbors. The tunneling constant J, which quantifies 
the microscopic tunneling rate between adjacent sites, is 
given by 



J ^ 



1 



i ( Vipn ■ V^Sfj+i ) -I- >pnV<fn+i 



. (35) 

with iffi being the wave function localized at site n. The 
on-site interaction as given by 



U = 



(36) 



is a measure of the interaction between atoms at one 
lattice site. The ground state of this Hamiltonian is a 
constant wave function V-'n = 1- Its excitation energy is 
given by e(q) = 2| sin {qn) \\/2JU , which yield the sound 
speed via Eq.lTI 



9^0 q 



(37) 



This result is consistent with the sound speed definition 
in terms of compressibility k and effective mass m* since 



1 



U : 



(38) 



In the following, we try to express J and U in terms 
of c and v. For a state well-localized at each lattice site, 
we can regard it as the ground state of the lattice well 
and describe the localized state (pfi with the ground state 
wave function of a harmonic oscillator. This approxima- 
tion immediately leads to an estimate of U . We obtain 



f/ = c (47r 



(39) 



AsC/; 



this indicates that the compressibility in the 



tight-binding limit is very similar to the weak potential 
limit: it depends on the dimensionality of the lattice and 
decrease with u in a non-exponential form. 

Mathematically, the time-independent Schrodinger 
equation for an atom in the cosine potential is a Mathieu 
equation. The theory of the Mathieu equation allows us 
to estimate J, which is given byjs^, HH 



(40) 



This result is very different from the weak potential limit: 
the effective mass ^ = -^^tj increases exponentially 
with ^/v. As a result, J should dominate the behav- 
ior of the speed sound. Combining Eq. ([39|) and Eq. (j40|) . 
we arrive at 

V, = 25/27r3/4cl/2 (4^2^)^/8^3/8 g^p(_4^)^ (4^) 



which shows the speed of sound decreases monotonically 
with V in an exponential form in all three dimensions. 
The sound speed in the tight-binding limit has a weak 
dependence on the dimension D of the lattice as D only 
appears in the prefactor of the exponential. 



IV. NUMERICAL RESULTS 

We have so far studied analytically the sound speed 
of a BEC in an optical lattice. In this section, we study 
the sound speed with numerical methods. Our numerical 
method allows us to find the sound speed for the whole 
range of lattice strength, particularly, intermediate lat- 
tice strength for which no apparent analytical approach 
can be found. 



A. Numerical methods 

As discussed in Section II, to compute the sound speed, 
one has to first find the ground state of the BEC system 
or the Bloch states in the vicinity of fc = 0. To find these 
states numerically, we expand the Bloch states in Fourier 
series 



N 

E 

m.n,l — — N 



i{mx-\-ny-\-lz) 



(42) 



where N is the cut-off. We find that TV = 5 is good 
enough for all dimensions. The Fourier coefficients 
{am,n,i} satisfy the normalization condition 



iV 

E I 

m,n,l — — N 



(43) 



Note that the Fourier coefficients {am,n,i} can be cho- 
sen as real. This fact greatly reduces the computation 
burden. 

The Bloch waves can be numerically obtained by vary- 
ing {am,n,i} so that the wave function (f)^ minimizes the 
system energy of Eq.®; the accuracy is checked by sub- 
stituting the solutions into the Gross- Pitaevskii equa- 
tion We use the standard minimization routine 
of MATLAB. The accuracy of the numerical solutions 
can be checked by substituting the numerical solutions 
into the time- independent Gross-Pitaevsku Eq.dU. Once 
the Bloch states <p^{r) have been obtained, we can com- 
pute the sound speed in two different methods. In one 
method, we calculate the Bogoliubov excitations e{q) of 
the ground state (poif) and find the sound speed of the 
BEC through Eq. (fT6|) . In the other method, we can 
calculate the effective mass m* and compressibility k, re- 
spectively, with Eqs.® and p8|) . Then the sound speed 
can be computed via Eq. (fT7|) . We have calculated the 
sound speeds with both methods and the agreement is 
excellent as expected. 
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B. Results and discussion 

We have computed numerically the sound speeds for 
all three dimensions for a wide range of lattice strength 
V and inter-atomic interaction c. The results are plotted 
in Figs lll2l & [3l respectively. Fig. [1] displays the sound 
speed in the one-dimensional case, which falls monoton- 
ically with increasing lattice strength. This is in agree- 
ment of previous studies [l6l. [TtI. llSj 
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FIG. 1: Sound speed for a BEG in a ID optical lattice via the 
strength of the optical lattice, (i) numerical result (solid line 
curve); (ii) analytical results for weak potentials (★ curve); 
(iii) analytical tight-binding results(B curve). 
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FIG. 2: Sound speed for a BEG in a 2D optical lattice via the 
strength of the optical lattice, (i) numerical result (solid line 
curve); (ii) analytical results for weak potentials (* curve); 
(iii) analytical tight-binding results(B curve). 

The results are different in two and three dimensions. 
As shown in Figs. [2]and[3l the relationship between the 
sound speed and the lattice strength depends crucially 
on the strength of interatomic interaction. In the two- 
dimensional case, when the interaction is above the crit- 
ical value, i.e. c > j, the speed of sound also decreases 
monotonically with increasing lattice strength (Fig[2I^d)). 
However, when the interaction is weak, i.e. c < 4, as 





(b) 


c=0.3| 










2.0 0.0 0.5 1.0 1.5 2.0 2.5 3.0 



FIG. 3; Sound speed for a BEG in a 3D optical lattice via 
the strength of the optical lattice, (i) numerical result (solid 
line curve); (ii) analytical results for weak potentials (★ curve); 
(iii) analytical tight-binding results(B curve). 



shown in Fig. [2Ia,b), a sound speed reaches a maximum 
at an intermediate strength of optical lattice. Fig[2fc) 
shows the transition point between the above two dif- 
ferent behaviors, where the sound speed changes almost 
does not change with weak lattice potentials. 

In three dimensions, the behavior becomes even richer. 
There exists a critical value of the interatomic interac- 
tion, c = i. When the interaction is smaller than this 
critical value c < ^, the sound speed first increases then 
decreases as the lattice strength increases (Fig. |3Ia,b,c)). 
This is similar to the two-dimensional case and was first 
noticed by Boers et al.^^. However, when the interac- 
tion is strong enough, i.e. c > ^, a new pattern is found. 
As shown in Fig. ^d), the sound speed can even oscillate 
with the lattice strength. According to our numerical re- 
sults, the oscillating behavior of the sound speed does 
not disappear until the interatomic interaction c reaches 
1 (c=l). 

Our numerical results are compared to our analytical 
results. As seen in Figs. [1] [2] and [3l our numerical re- 
sults agree very well with our analytical results (★ curves) 
in the regime of weak potentials. For strong lattices, our 
numerical results also agree well with the tight-binding 
results (■ curve in Figs. [H Ula) and[3l^a)) for weak in- 
teractions. However, for strong interaction in two and 
three dimensions, there exist discrepancy between the 
tight-binding results and our numerical results. This dis- 
crepancy is very large especially for three dimensions as 
shown in Fig. [31[b,c,d). The agreement can only be best 
described as qualitative. 

This big mismatch is nevertheless expected as both 
Eg. ([55]) and Eq. (|i(I)) are derived without the considera- 
tion of interaction between atoms. It can be explained 
by noticing that the interaction can strongly modify the 
ground state wavefunction localized in each well and 
greatly enhance the tunneling between the adjacent wells 
in two and three dimensions. As in Eq. (j38p . the tun- 
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neling rate J is related to the effective mass m*, which 
is very sensible to the behaviour of the wave function 
in the region of the barriers. In Figs. IHal)-(a2), Figs. 
[SIal)-(a4) and Figs. [nilal)-(a4), we show the effect of 
the interaction c on the effective mass m*. The effect is 
different in different dimensions. For one dimension, as 
shown in Figs. [Hal)-(a2), the effective mass m* increases 
greatly with the lattice strength for both weak and strong 
interaction. In 2D, the interaction has much stronger in- 
fluence on TO*. If we compare FigEl^al) and Figl5I^a4), 
the influence is order of magnitude different. In 3D, the 
interaction c affect the effective mass to* most strongly. 
As shown in Fig. [H^al), for weak interaction c — 0.01, 
the effective mass increases to ^ 150 aX v — 2. In Fig. 
[6l^a4) for c = 0.7, the effective mass m* is only ~ 3 at 

V — 2. This order of magnitude difference explains why 
there are large disagreement between our numerical re- 
sults and the tight-binding model for the speed of sound 
in FigEl^b-d). 

The dependence of sound speeds on the lattice strength 

V is largely expected from our analytical results for the 
two limiting cases of weak and strong lattices. We have 
shown that in the weak lattice limit, the speed of sound 
can either increase or decrease with the lattice strength 
while in the strong lattice limit the speed of sound always 
decreases with increasing lattice strength. Naively, one 
would expect that the speed of sound either decreases 
monotonically with the lattice strength v or develops a 
maximum at certain intermediate value of v. This is 
exactly what we have seen in Figs. [TJ [5] and [3] except in 
FiglJJd) where we see two local maxima. 





0.0 0.5 1,0 1.5 2,0 0,0 0,5 1,0 1,5 2,0 



FIG. 4: Inverse compressibility kT^ , effective mass m* and 
on-site interaction energy U for a BEG in a ID optical lattice 
via the strength of the optical lattice. 

To better understand the behavior of the sound speed 
as a function of the lattice strength u, we have also com- 
pute numerically the effective mass m* and the com- 
pressibility K and the results are plotted in Figs l4l5i &l6l 
It is clear from these figures that the compressibility n 
has different dependence on v in different dimensions. 
This agrees with our analytical results in the last sec- 




0,0 0,5 1,0 1,5 2,0 0.0 0.5 1.0 1.5 Z.O 



V V 



FIG. 5: Inverse compressibility effective mass m* and 

on-site interaction energy U for a BEG in a 2D optical lattice 
via the strength of the optical lattice. 




V V 



FIG. 6: Inverse compressibility , effective mass m* and 
on-site interaction energy U for a BEG in a 3D optical lattice 
via the strength of the optical lattice. 



tion. However, we notice that the increasing rate of k -"^ 
with V is quite close in all dimensions. 

The situation is different for the effective mass m*. In 
the last section, we have shown that the effective mass 
TO* has the same dependence on v (see Eq. (|5T|) ) in all 
dimensions. However, it is true only in the limit of weak 
lattices. As seen in the right columns of Figs l4l51 
the effective mass to*, as a function of behaves very 
differently in different dimensions. In particular, the in- 
creasing rate of m* with v in one dimension is orders of 
magnitude larger than the increasing rate in three dimen- 
sions. The two-dimensional case is right in the between. 
Since the sound speed is the result of competition of m* 
and K, the relatively small increasing rate of to* with v 
allows the sound speed oscillates with v in 3D. 
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V. EXPERIMENTS 

The speed of sound of a BEC in an optical lattice may 
be measured with the similar technique that was used in 
Ref. ^13*1 to measure the sound speed of a BEC in a trap. 
Some complication is expected due to the periodic modu- 
lation of the BEC density. Another possible method is to 
employ Bragg spectroscopy [H, HI, [s^l to the excitation 
spectrum. The speed of sound can be extracted from the 
slope of the linear part of the excitation spectrum. 

In typical experiments to date, the relevant parame- 
ters are as follows: for a BEC in three-dimensional op- 
tical lattice[ll[, the atom occupancy per lattice is of the 
order of (n,) « 1 ~ 3, no = 1.3 - 3.9 x lO^^m-^, 
fci = 2ti/\l = 7.37 X 10^TO"\ and as = 5Anm; 
for a BEC in quasi- two-dimensional optical lattice [s^. 
the atom occupancy per lattice can reach (rii) ^ 170, 
no = 3.6 X 1020to-3, fc^ = 2n/XL = 7.37 x 10Sn-\ 
and flj = 5.4nm; for a BEC in quasi-one dimensional 
optical latticefH, (n,) w 1000, no = 2.8 x lO^^rn^^, 
fci = 2tt/Xl = 7.9 X lO^m^^, and = 5.4nm. These 
parameters correspond to c = 0.08 for ID optical lattice, 
c = 0.11 for 2D optical lattice, and c = 0.004 - 0.012 for 
3D optical lattice. The depth of optical lattice Vq can be 
changed from OEn to 32£'i^'Tl^; it means that our v can 
be changed from and 2. To our knowledge, the highest 
atomic density without lattice is no = 3 x lO^^m"^ for 
sodium[36]. For this high density, we have c = 0.22 with 
ki = 1.07 X 10^ m^^ 4]. However, this is rather idealis- 
tic. The other possible way to increase c is to tune the 
scattering length Og with Feshbach resonance [stI [ssj. 



VI. CONCLUSIONS 

We have studied the speed of sound, compressibility, 
and effective mass of a Bose-Einstein condensate in an 
optical lattice both analytically and numerically. Special 
attentions have been paid to the effect of the depth of 
the optical lattice v, the interatomic interaction c and 
the dimensionality D on the sound speed. Our investiga- 
tion shows that the sound speed depends strongly on the 
strength of the lattice. In the one-dimensional case, the 
speed of sound falls monotonically with increasing lat- 
tice strength. The dependence becomes much richer in 
two and three dimensions. In the two-dimensional case, 
when the interaction is weak, the sound speed first in- 
creases then decreases as the lattice strength increases. 
For the three-dimensional case, the sound speed can even 
oscillate with the lattice strength. These rich behaviors 
can be understood in terms of competition between com- 
pressibility and effective mass. Our analytical results at 
the limit of weak lattices also offer an interesting perspec- 
tive to the understanding: they show the lattice compo- 
nent perpendicular to the sound propagation decreases 
the sound speed while the lattice components parallel to 
the propagation increases the sound speed. 
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APPENDIX A: PRELIMINARY NOTATIONS 

Suppose / (r) to be a periodic function with the peri- 
odicity of R, given by 

/(r) = /(f+i?), (Al) 

with 

R = midi + 11120,2 + m-afls, (A2) 

where f is the position vector, oi, a2, and are any 
three vectors not all in the same plane, and toi, 1112, and 
rna ranges through all integral values. Corresponding to 
ai's, there exist a set of reciprocal vectors fej's such that 

flj • bj = 2TTSij. (-A-3) 

We can expand the periodic function / (r) as its Fourier 
coefficients ^nif ) as defined by, 

/ (f) = ^ (/) exp {in ■ f } , (A4) 

n 

with 

^n{y)^\z f drf (r) exp {-in-r}, (A5) 

" J cell 

and 

n = nibi + n2b2 + n^bs. (A6) 

where the nj are integers. In the integration, Q is the 
volume of the primitive cell and the integration is over a 
primitive cell. 

APPENDIX B: SOLUTIONS OF THE 
GROSS-PITAEVSKII EQUATION IN THE WEAK 
POTENTIAL LIMIT 

The time-independent Gross-Pitaevskii (CP) equation 
in the three-dimensional case can be written as 

-^v'^(?') + c|V'(^lV(r) + l^a,.(r)^(f) =^V(r), (Bl) 

where Var{r) is the periodic potential with the periodicity 

of R, 

Var - Var {f + R) ■ (B2) 
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The Bloch-wave solutions of the GP equation (IB1[) reads 

^Pir) = (/.,r(f)e*n (B3) 

where k is the Bloch wavenumber and 4'j:{r) is a periodic 
function with the same periodicity of Eq. (jB2p . Substi- 
tuting Eq. I|B3|) into Eq. (|Bip . we have the following 
equation for each Bloch wave state (t>^{r) 

-\{V + i^fh{^ + c|(/.^:(r)|'0,-(f) 

+K.(r)0,-(f) = /i (fc) (B4) 

The set of eigenvalues then forms Bloch bands. 

Besides the GP equation (|Bip , the Bloch wave function 
is also subject to the normalization condition given by 

i / dm'' = 1, (B5) 

" J cell 

which is equivalent to 

^g(|</.|') = 1. (B6) 

For convenience, we have dropped the suffix k and the 
coordinate vector r in (jij: (r) 

Expanding cj) in terms of the potential strength as 

= 0W+0(i)+0(2)+... , (B7) 

we get the zeroth, first, and second order forms of Eq. 
()B6| . respectively, 

n 

n 

+^a(</)W)^i(0W)) =0, (B9) 

n 

There is still an arbitrary phase in the above wave func- 
tions, which satisfy both the GP equation and the nor- 
malization condition. Therefore, we may impose a third 
condition 

^ J e R, (Bll) 

So that the Bloch states can be uniquely determined. 

Before solving the GP equation (|B1[) . we have to set 
forth another two specifications. First, we are only con- 
cerned with Bloch states at A: = 0. In this case, we rewrite 
Eq. (|B4)) as follows 

-lv^<|) + c\(|>U + Var{r^<|> = ^^<|>, (B12) 



where we dropped the suffix and the coordinate vector 
r in for convenience. Expanding and fi in terms 

of the potential strength, we get the zeroth, first and 
second order forms of Eq. (jB12[) . respectively, 

_iy2^(0)+c|<^(")|V°) (B13) 

_ iv20(i) + c (2 1 (/.(o) I '0(1) + 

+K.(r)0(") = + (B14) 

+20(O)|</)«|' + 0(O)VW') +K,(r>« 

= /i^°V^'^+M^'V^'^+M^'V^''^. (B15) 

Second, we are only concerned with the cases in which 
the external potential Var{r) is symmetric in each cell, 
or in other words, Var{r} is an even function. Combining 
it with the condition that Varif) is a real function, we 
immediately have 

■^n{V) = .^-niV) e R. (B16) 

In the following, wc will solve the GP equation for ob- 
taining the normalized Bloch state a,t k — 0. 

1. The zeroth-order correction of the GP equation 

From Eq. (jB13[) , we get the zeroth-order wave function 
and chemical potential, respectively, 

<^(0) = 1, ;i(")=c, (B17) 

which automatically meet the normalization condition 
(|B8l) . 

2. The first-order correction of the GP equation 

Substituting Eq. (|B17p into Eq. ^9^, we have 

=^o('^^'^) + ^g(0^'^)-O. (B18) 

From the phase condition (jBlip . we know that ^f^{(j)'^^^) 
is a real number, and therefore 

^g(0(i))=O. (B19) 

Substituting Eq. (|B17p into Eq. (|B14p . we have 

+^n{V)=fi^'^6^s. (B20) 
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Plugging Eq. (|BT9| into Eq. (|B20)) and letting n = 0, we 
get the first-order correction of the chemical potential 

Ai(i) = ^^(V). (B21) 

Taking complex conjugates on both sides of Eq. (|B20p 
and replacing —ft with n, we obtain 

+,^n{V) - A*(i)<5,,ff. (B22) 

The unique solution of Eqs. (|B20p and (|B22p in the case 
of n 7^ reads 

^,(</>W)=^,(0W*) = -^^?gL, n^O. (B23) 

^\n\'^ + 2c 

From Eqs. (|B19p and (|B23|1 . we know 

^«(</)W)=^_,-j(^(i))eR, (B24) 
which means that 4>^^\r) is a real even function. 



Since the correction to the system energy is second order 
in the potential strength, it is sufficient to consider the 
first-order correction of the Bloch state, there is no need 
of calculating the second-order correction of the Bloch 
state. Based on Eq. (jB29[) . we can easily obtain the 

system energy E ^fcj near A; = 0, up to the second-order 

correction. 



APPENDIX C: ANALYTICAL EXPRESSION OF 
SOUND SPEED BASED ON EQ. (17) IN WEAK 
POTENTIAL LIMIT 

The aim of this section is to calculate the compress- 
ibility K and the effective mass m* as a function of the 
interatomic interaction c and of the depth of the arbi- 
trary periodic potential Var {r)- Using these quantities, 
we will calculate the velocity of sound. 



3. The second-order correction of the GP equation 

Plugging Eqs. (|BT71) and (|BT91) into Eq. (|BTo1) . we 
obtain 



1. Compressibility k and effective mass m* 



Plugging Eqs. (IBTTI) . (IB2T]) . and (lB28l) into Eq. P^. 
we obtain the analytical expression of compressibility k 
=Fo(</>(2^) + ,^Q{(t>^^^)* + ^ .^ni^)^^^)^ = 0. (B25) in the weak potential limit. 



Plugging Eqs. (|BT71) and (|Bl9)) into Eq. (|BT51) . we obtain 



1-E 



16|n|^ 



v 



n#0 (^\n\' 



4c 



(Cl) 



In the case of n = 0, we have 

= A.W^g(</.«)+/.(2)<S^5. (B27) 



To calculate the sound speed, we also need calculating 
(B26) the effective mass m*. Substituting Eq. I|B30|) into Eq. 
([S]), we obtain the analytical expression of effective mass 
along a given direction indicated by a unit vector f. 



16 |n • 



(C2) 



Plugging Eqs. im]), (EMI and iB25| into Eq. jB27|, ^ ^e^^f « j^f * effective mass along each axis x, 



we 



obtain the second-order correction of the chemical ^' J z 



potential 



n^^O V 2 I 



2c) 



(B28) 



and 



To complete the calculation of the sound speed, we 
also need to calculate the system energy near k = 0. 
This can be obtained in terms of the effective potential 
c|(/)p -|- Var — M seen by each atom. We view our system 
as a noninteracting gas in the effective potential. 



Veff (0 



n H -I- 4c 



Vir). 



(B29) 



1 



1 

m* 



1 

TO* 



1-E 

n#0 



1-E 



16 |n • afj^ 



1-E 









|n| 




+ 4c)^ 




16 |n- 










|n| 




+ 4c) 




16 |n • 


^,2 
^1 



n^o \n\ \n\ +4:C 



(C3) 



(C4) 



(C5) 
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Plugging Eqs. ([02]) and ([CT]) into Eq. dTT]), we arrive 
at the analytical expression of sound speed labeled by 
along a given direction f, 



(4c+ |n|2)3 



E 



In • 



|n|2(4c+ |n|2): 



r}^B2(T/). (C6) 



Plugging Eqs. jell, del, Jen, and ^ into Eq. 
^T7\ . we also obtain the analytical expressions of sound 
speed along each axis x, y, and z, labeled by v^^, Vgy, 
and Vsz, 



Vsx = VC 



l^„(4c+|nPf 



^JnP(4c+|npy 



and 



(4c+ |ri|2)3 



E 



and 



8V^{E 



(4c+ |n|2)3 



E 



In • 



|n|2(4c+|n|2): 



:]^n\V). (C9) 



In the following, we consider a special case, i.e. ai, 
02, and 03 are chosen along x, y, and z respectively; we 
also suppose |ai| — |a2| = lasj — 27r, without loss of 
generality. In this case, the sound speed of Eqs. (|C7p . 
(|C8|) . and (|C9p can be simplified into, 



(4c + n2)3 



i:_^,i^^.Wiy)^ (CIO) 



and 



8Vs{E 



(4c + n2)3 



E 

ft 7^0 



n2(4c ■ 



}^^iV), (Cll) 



I- ^ (4c + n2)3 



E 



:}:^n^{V). (C12) 



APPENDIX D: ANALYTICAL EXPRESSION OF 
SOUND SPEED BASED ON EQ. (16) IN WEAK 
POTENTIAL LIMIT 

As shown in Section II. B, there are two equivalent ways 
to calculate the velocity of sound; one is based on Eq. 
(16), the other comes from Eq. (17). The aim of this 
section is to calculate the analytical expression of sound 
speed based on Eq. (16) from the another angle, by di- 
rectly solving excitation energy e{q). 



1. The matrices P, Q, S and T 

According to the Bogoliubov theory, the excitation en- 
ergy e{q) of the EEC in the Bloch state at fc = can be 
obtained by solving the following eigenvalue problem 



with 



M 



and 



^{k + q) c(j)j:^ 
ccpf ^{~k + q) 



1 

-1 / ' 



(Dl) 



(D2) 



(D3) 



where ^{p) is defined as 

^{q) = -^(V + + Kr(r) - + 2c|0,-|'. (D4) 

By a similarity transformation, we can transform CTzM 
into a numerical matrix P without changing the eigen- 
values. The new matrix P can be represented in a block 
form 



(D5) 



where each block is actually a 2 x 2 matrix and 
rfi and n take values ranging from (—00, —00, —00) to 
(-I-00, -l-oo, -l-oo). For convenience, we abbreviate the di- 
agonal blocks Tf^ft in Eq. (jD5[) as Sf^^ and consequently 
T^fi denotes solely those non-diagonal (to 7^ n) blocks. 

For we are only concerned with the case of A: = 0, CTzM 
can be simplified into 



(D6) 
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In this case, we have 



with an, bfi, c^n and determined by 



(D7) 



an = + q}^ + ^o{V) - + 2c^o(l0l'), (D8) 



(D9) 
(DIO) 

dmfl =C^m-n(</'^). (Dll) 



bn = C^o 
Cfhn — ^ m— n(^) ~t~ 2Ct^ n(|0| ): 



We define the matrix Q as follows 
Q = P - £/, 



o. The zeroth-order approximation of e 

The zeroth-order form of Eq. (|D13[) is 

(det Q)<"^ = det Q(°) = det(p(°) - = 0. (D23) 

From Eqs. (jPlSp . we know that all the T^l are zero ma- 
trices, and therefore the matrix is block diagonal as 
represented in Eq. (jD5|) . Consequently, the eigenvalues 

of P^"-* are the collection of the eigenvalues of each s'^\ 
The zeroth-order approximation of e is hence the positive 

eigenvalue of Sq 



(0) 



6f 



= 0. 



(D24) 



(D12) rpj^g positive solution of Eq. ([DM]) reads 



where e represents the lowest elementary excitation. 
Consequently, we have 



.(0) 



4 



(D25) 



det g = 0. 



(D13) 



Now, we compute the elements a,-j, bfi, and d^^j from 
Eqs. jDSl - (IdTT|) as follows 



a:-' - l{n + q)^ + c, bf ^c, (D14) 



= 0, = 0, 
(1) _ \{m - nf - 2c 
^(m~n) + 2c 



(D15) 
(D16) 

^,n-niV), (D17) 



2c 



mn I 



,^ ^,2 . ^rn-niV), (D18) 

2 (to — nj + Ic 



.(2) 



= 0. 



1^2 



r^«(l^)^ (D19) 



(D20) 



With the value of e^*^^ , we calculate the determinant of 
each diagonal block of the matrix Q^"^ 

det(4°) - e(«)/) = {^-f + c)' - Q(n + + c)', 

(D26) 

(D27) 



which are denoted by for convenience 



This result will be useful in the following sections. 



6. The first- order correction of e 

We can conclude that the first-order correction of e 
vanishes as 



(D28) 



For all these quantities are real, we may drop the asterisks 
* in Eq. (jPTl) . 



c. The second-order correction of e 



The second-order form of Eq. (|D13p reads 



2. The lowest elementary excitation in the weak 
potential limit 

We expand the matrix Q and e in terms of v as follows 



e = e(o)+e(i)+£(2) + 



(D21) 
(D22) 



The aim of this subsection is to calculate e*^°-', e^^-* and 
e^^^ by expanding Eq. (ID13|) into its zeroth, first and 
second order form. 



V (^Qij ) 

"J 



2 ^ \dQijdQki 

ijkl ■' 



Here we introduce the "second cofactor matrix" of Q, 
which is denoted by Q and whose elements are defined 



as 



(D30) 
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By this notation, we reduce Eq. (|D29p into 



E 



2 E QijMlQij^Qk ^ 
ijkl 



0. 



(D31) 



We start by computing the first term on the left hand 
side of Eq. (jD31|l . Keeping only the non- vanishing terms 
and noting that = 0, we have 



(0)£(2) 



r, (0) (2) 



n 



(2e 



(0)5(2) 



rh^n 
m/0 



Then we proceed to compute the second term on the 
left hand side of Eq. (|D3ip . It should be noted that any 
Qij is one of the elements of the matrix U^fi which is 
defined as 



b,n 




drhfi 












bn 




-bn 


-an - 



(D33) 



The second term on the left hand side of Eq. (jD3ip 
could be computed in a routine way by computing the 
second cofactor matrix of U^nn- However, we have found 
a much more convenient method to compute this term 
which is shown as follows. Since we have 



r(l) - 



-o^-e^''=b^=-b^-0, (D34) 



the only non- vanishing Q^j^ are those in T^l, or more 
specifically, c^l, -c^l, d^l and -d^l. We further as- 

^(0) 

sert that in order to get a non- vanishing Q^j , the cor- 
responding Qij and Qki must be in the same U^n- This 
assertion can be confirmed by some routine proof which 
we shall not elaborate here. From the above assertion, 
we obtain the following expression 



E 

ijkl 



m^n 4 
m,n ijkl—1 



{U^nf^l,{U^n)[]\u^nt^ Y[ %.(D35) 



Keeping only the second-order correction terms, we have 

4 

ijkl — 1 



mn nm 
,(0)1,(0) Jl) ,(1) 



m n mn nm 
,(0)i.(0)Jl) j(l) 



m n mn nm 



n m nm mn m n nm mn 



n m mn nm 

-2a^?a^^U^llS}l 

m n mn nm 



m n mn nm 

2bWd^lU}l, 

m n mn nm^ 



which is denoted by '^^n for convenience. Since both 

mn 

for 



and d^^^ are symmetric in rh and n, the same is true 



mn ^'^nm- 



(D37) 

Plugging Eq. p36p into Eq. p35p and noting that 
= 0, we obtain 

E - 2 ;^ 5^ n (D38) 

Plugging Eqs. (|D32l) and (|D38|l into Eq. (|D3T|) . we have 
(2s(°)s(^) - 2afaf) R + E 11 ^'^ = 

(D39) 

from which we finally arrive at the value of e^^-* 



1 

n^o V V2I" 



2 + 2c)' 2^^^ 



(D40) 

3. The speed of sound in the weak potential limit 

The speed of sound along any given direction f in a 
BEC is defined as 



Vn = |V,-e(9)|-_^g+ 



(D41) 



Let's consider a special case, i.e. oi, 0,2, and 03 are cho- 
sen along each axis of reference system, x, y, and z, re- 
spectively; and without loss of generality, we suppose the 
periodicity of the periodic potential along each axis to be 
2tt. We are particularly interested in the speed of sound 
along one of the three axes, for example x-axis 



de 



(D42) 



Plugging Eqs. ((D25)) . (|D28|) and (lD40l) into Eq. ((D42|) . 
we finally obtain the analytical expression of the sound 
speed along each axis, 

Vsx = VC+ 2^ — — -—3 [V), 

(D43) 



n^o 



and 



Jsy 



n^o 



|n|2(4c-f |n|2)3 



and 



(D44) 



v^+E 



8\/c [(n? + ni)|np - 4cn§ 
|nP(4c-h |n|2)^ 



(D36) 



(D45) 

It can be easily proved that Eqs. (|D43P - (|D45|) can be 
deduced into Eqs. ((CT0| - ((CT2|) . 
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APPENDIX E: A SPECIAL EXAMPLE 



Eqs. (7) and (8), respectively, 



Suppose that the arbitrary potential of Var (f) is chosen 
to the special form of Eq. ^ 

Viatt {r) = V (cos X + cos y + cos z) . (El) 

In this case, there are only six non-vanishing Fourier co- 
efficients 

■^fi{V) = - (^rl, (±1,0,0) + <5»i,(0,±l,0) + '^n,(0,0,±l)) ■ (E2) 



Submitting Eq. (|E2| into Eq. p43|l . we have 

y,^^c(l + p^^:^vA , forD = 3. (E3) 
V (4c +1) J 

With the similar calculations, we can also obtain the an- 
alytic expressions of sound speed in the optical lattice of 



v.,^. = V c 1 



16c 



(4c + iy 



:V^\jorD^l, (E4) 



'-^^vAjorD = 2. (E5) 



Combining Eqs. (jE3p . (jE5|) and (|E5|) together, we arrive 



= Vc 1 



4(L> - 1 - 4c) 
(4c +1)^ 



(E6) 
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